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We experimentally examine the dynamics of two-particle collisions occuring on a surface. We find 
that in two-particle collisions a standard coefficient of restitution model may not capture crucial 
dynamics of this system. Instead, for a typical collision, the particles involved slide relative to 
the substrate for a substantial time following the collision; during this time they experience very 
high frictional forces. The frictional forces lead to energy losses that are larger than the losses 
, due to particle inelasticity. In addition, momentum can be transfered to the substrate, so that 

q_) ■ the momentum of the two particles is not necessarily conserved. Finally, we measure the angular 

Fj^ | momenta of particles immediately following the collision, and find that angular momentum can be 

. lost to the substrate following the collision as well. 
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I. INTRODUCTION 



Dry granular systems have generated much interest recently in the physics and engineering communities, both 
for fundamental understanding and for direct applications These systems are important both in nature (e.g. 

£H , avalanches) and in industry (e.g. pharmeceuticals and grain elevators). 

Particles in such systems are typically considered to interact only through interparticle collisions, i.e. repulsive 
O , contact forces. Experiments that can yield quantitative data for velocities, collision rates and other useful quantities 
are often performed in two dimensions. In order to allow reasonable motion of particles in such an experiment, the 
particles must either be free to roll or they must be levitated, for instance by air flow. Here, we consider particles 
rolling on a smooth flat surface. We note that there are then two types of friction that the particles experience when 
in motion. The first, rolling friction, occurs when the particle is moving without sliding on the substrate; its effect 
is relatively weak, with a coefficient of friction on the order of 10~ 3 Q|. Rolling friction affects individual particles 
independently of collisions; it tends to damp motion slowly over time. It also affects the mobility of particles on 
the surface. For example, segregation occurs when particles of differing surface properties are shaken on a smooth 
surface The second type of friction affecting particles is sliding friction. This occurs when the contact point 
of the particle and the surface is not instantaneously at rest. Sliding friction can occur when particles undergoing 
collisions experience frictional frustration, i.e. when it is impossible to maintain nonsliding contacts between colliding 
particles and the substrate. During a collision the contact force between the particles is much greater than the force 
of gravity, so some sliding on the substrate will occur. Sliding friction is much more dissipative than rolling friction, 
with a coefficient of friction on the order of 10 _1 In the experiments described here, the sliding interactions with 
the substrate are the predominant mechanism for energy loss. The sliding of particles following a collision leads to 
an energy dissipation rate that is ~ 10 2 times greater than dissipation from rolling friction. The time over which 
particles slide is typically relatively long, ~ 0.05 — 0.1s. Hence, the effective time over which a collision influences the 
dynamics of a particle is much longer than the actual contact time of ~ 10~ 5 s 0. After a pair of particles has stopped 
sliding, the momentum of their center of mass (in the lab frame) need not be the same as the before-collision value. 
These features have tremendous importance on dynamics of systems rolling on surfaces, but have been relatively little 
explored experimentally. However, recent related theoretical and numerical work has been conducted by Kondic ||. 

The purpose of this paper is to examine in detail some important aspects of the dynamics of two-particle collisions 
which occur on a surface. We begin by briefly describing the measurement apparatus used to follow the particles' 
motion, and then discuss the surface effects. 



II. APPARATUS 



The particles used were 2.38mm steel balls, which moved on a flat aluminum surface. The aluminum was black 
anodized to improve visual contrast between the steel spheres and the background. The apparatus was illuminated 
from nearly directly overhead; with this lighting, each metal sphere produced a single bright spot near its highest 
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point due to the reflection of the overhead light. In order to track the centers of individual spheres over time, we used 
high speed video at rates of 250 frames per second. We then used particle tracking techniques to follow the particles. 
We began by finding the positions of the centers of all particles within a video frame, identified by the brightest points 
in the image (the local maxima in the brightness field). Although there were some secondary reflections between 
neighboring balls, these reflections were much less bright than the primary reflections, and they could be eliminated 
easily. By following the positions of individual particles from frame to frame, we obtained trajectories, velocities, and 
other time- varying quantities of interest. 



III. PARTICLE-SUBSTRATE DYNAMICS 



A. Rolling friction 

The simplest effect of motion on a substrate is rolling friction, and we consider this effect first. The frictional force 
from a single sphere rolling on a substrate is usually modeled by 

Ff r — ^Fn, (1) 

where Fn is the normal force at the sphere-substrate contact and /i r is the coefficient of rolling friction. 

We have carried out measurements of the frictional force on a single steel sphere rolling on the aluminum substrate 
described above. The sphere was tracked as described in Section We determined its acceleration by dividing the 
change in velocities between two frames by the time between the frames; the resulting acceleration versus velocity is 
shown in Fig. [|. The solid line in the figure corresponds to a least-squares linear fit to the data. We see that the rolling 
friction for this system is velocity dependent, with higher frictional force at higher velocity. This tends, in principle, 
to make velocities in rolling granular systems become more uniform. To a reasonable approximation, the acceleration 
due to rolling friction which a particle experiences is a — —Cv — D, where C = 0.135s" 1 and D = 1.03cm/s 2 . Using 
a = ji r g, with a typical acceleration of a = — 2.5cm/s 2 , we find that fi r ~ 2.5 • 10~ 3 , which is comparable to that 
reported by Kudrolli et al. 0] for steel balls rolling on a Delrin surface. 



B. Sliding friction 

While this rolling friction has a dissipative effect over long times, there is another, stronger, mechanism for energy 
loss to the substrate: sliding friction. We find that sliding friction with the substrate immediately after a collision 
plays a particularly important role in the system dynamics. In order to investigate this effect, we consider the collision 
of two particles on a substrate. We first review the textbook example of two particles colliding in free space, which 
we assume is two-dimensional, and then compare this to experimental observations when the motion occurs on a 
substrate. 

In the standard case of two inelastic frictionless particles colliding in free space, i.e. with no substrate, the collision 
is described by conservation of momentum and by an energy loss given through the coefficient of restitution, r. We 
introduce the following notation to describe this process. The initial momenta of the two particles are given by pu 
and p2i, the final momenta by p\f and P2f, and Pi = pu + pu- The direction of the vector connecting the centers of 
mass of the two particles at the time of the collision is h. The relative velocity of the particles in the n direction after 
colliding is a fraction r of their initial relative velocity, while the relative velocity tangential to n is a fraction s of its 
initial value. Thus, 

Pif + Piif = Pu + P2i, (2) 
Plfn - P2fn = -r(piin - P2in), (3) 
Pi ft ~ P2ft = s(p U t ~ P2it), (4) 

where the subscripts n and t refer to the directions parallel and perpendicular to n, respectively. We take s = 1, the 
simplest case describing an inelastic collision. 

This model is usually used in modeling granular systems However, it does not accurately reflect the dynamics 

of two rolling particles colliding on a surface. When two rolling particles collide, there are three contact points: each 
particle with the substrate, and the particles with each other. In general, these contact points are frictional. This leads 
to rotational frustration and, after the collision, to sliding fl. The following simple argument shows why the particles 
are likely to slide on the substrate following a collision. During a collision the frictional force between the particles 



2 



competes with the frictional forces between the particles and the substrate. If the static friction coefficients at all 
contacts are comparable, the frictional force will be greatest where the contact forces are greatest. The interparticle 
contact force is F p _ p ~ Ap/At, where Ap — mAv is the momentum change of a particle and At contac t ~ 10 _5 s || 
is the contact time for a hard-particle collision. The contact force for a particle with the substrate is F p - S = mg, 
and the ratio F p - p / F p _ s = Av/(g At contact)- Thus gAt conta ct ~ 10 -2 cm/s for hard metal spheres defines a crossover 
velocity, with sliding on the substrate occuring for Av > gAt con tact- If sliding has been initiated in the collision, 
a finite time and distance is required after the spheres separate before dynamic (sliding) friction slows the spheres' 
sliding motion. They will eventually reach a point where sliding stops and the particles are simply rolling. During 
this time, both the direction and the speed of the particles change significantly, as detailed below. 

We have investigated this effect experimentally by rolling one ball at an identical stationary ball, and by tracking 
their motion before and after the collision. Fig. ^ shows a typical set of trajectories in such a two-particle collision. 
The moving ball (in the lab frame) enters from the left of the image, hitting the stationary ball. Since it is difficult 
to produce a perfectly head-on collision, the incoming ball strikes the stationary ball slightly off-center. Immediately 
after the collision the two balls behave almost as though there were no surface interactions (see inset). Somewhat 
later, the particles begin to show the influence of the substrate as they change direction and speed. 

During the time between the collision and the time when the balls begin rolling without sliding, both the direction 
and the speed of the balls change due to sliding. This is shown in Fig. y, which gives the distance between the two 
balls shown in Fig. || over time. The collision occurs at t c w 0.06sec. For times t before and well after the collision, 
the separation s varies nearly linearly with time, indicating that the balls roll with nearly constant velocity in these 
periods. By contrast, during the ~ 0.06sec immediately following the collision, the interparticle separation varies 
nonlinearly in time. This indicates a regime in which the two particles experience dynamic, or sliding, friction with 
the substrate. We denote the time following the collision at t = t c and before the particles start rolling without sliding 
at t = t r as the "relaxation time," tr — t r — t c . We define t c as the time at which the particles' centers of mass are 
closest together, and t r as the point in time after the collision at which a particle begins moving with nearly constant 
velocity. r r was typically 0.05 — 0.1s in the systems we studied, which is very large compared to the time the particles 
are in contact, roughly 10~ 5 s ||. After a period of time equal to tr has elapsed, each particle has nearly constant 
velocity, affected only by rolling friction. 

From collision data we can determine the coefficient of restitution r, as defined in Eq. ||by examining the velocities 
immediately preceeding and after the collision, but before sliding friction has had significant effects. Fig. |] gives 
a histogram of data obtained for a number of measurements of r obtained this way. To produce these results, we 
measured velocities immediately before and within 0.01 seconds after the collision. We find an average value of the 
coefficient of restitution to be r = 0.85 ± 0.11, with no obvious dependence on the velocity of incoming particle. This 
value is similar to the value of r reported by others for steel-on-steel collisions: r = 0.93 in ref. [Q and r = 0.90 in 
ref. §. 

Immediately after the collision, the relative angle of the particles' new directions is also close to what one would 
expect for an elastic collision between two equal-sized spherical particles with one initially at rest. Fig. [| shows a 
typical example. In a collision between two identical spheres of radius R, with coefficient of restitution r and impact 
parameter b, the angle between the directions of the spheres' motion after the collision is given by 
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Then as r — ► 1, a perfectly elastic collision, cos8f = and Of = 7r/2, provided that b/2R ^> 1 — r. A typical value 
of the impact parameter in these experiments is b/2R ~ 0.25. The collision in Fig. occurs at time t « 0.06sec, 
indicated by the vertical dotted line; at this time, the relative angle between the particle velcities is near tt/2. As 
sliding friction begins to affect the particles and they are accelerated or deceler ated, the angle between the velocities 
decreases. The direction of the acceleration is discussed in detail below, in Sec. IIID| . 

Kondic has investigated a model for two particles colliding on a surface that includes both the interaction between 
the particles via a collision and the interaction of the particles with the substrate through friction || . For a system 
consisting of a moving particle hitting a stationary particle head-on, he predicts velocities and relaxation times of each 
particle after the collision. If the initial velocity of the moving particle is v , the final (i.e. purely rolling) velocities 
of the initially stationary and initially moving particles are V\ f and «2 / respectively, and the relaxation times of the 
initially stationary and initially moving particles are T\ and T2 respectively, then: 
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where r is the coefficient of restitution of the particles, \x k is the coefficient of kinetic friction of the particles with 
the substrate, g is the acceleration of gravity, and / is the moment of inertia of the particles (/ = 2mR 2 /5). C is a 
measure of the transfer of angular momentum between the particles during the collision, such that immediately after 
the collision 
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where u>i, for i > 0, is the angular velocity of each particle and oj is the angular velocity of the incoming particle 
before the collision. 

We fitted data for experimentally determined final velocities v\ and v 2 versus v a and relaxation times t\ and r 2 
versus v a to these predictions by minimizing the squared deviation of the model from observed data. Specifically, we 
minimized 
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with fitting parameters r, and C. Here a n represents the experimental uncertainty in the variable n. We found 
that r = 0.903 ± 0.008, fi k = 0.232 ± 0.023, and C = 0.347 ± 0.008 in our experiments. Note that this result for the 
coefficient of restitution r is consistent with, but much more precise than, the value of r = 0.85 ± 0.11 determined 
from Fig. [|. Results of these fits can be seen in Fig. § and Fig. @. 

We conclude that for a two-particle collision on a substrate, the picture of an instantaneous normal coefficient of 
friction is inaccurate, and may not be particularly useful. Without surface interactions, the relaxation times are r = 
and the final velocities are v\jv = (1 + r)/2 and v 2 /v = (1 — r)/2 (represented by the dotted lines in Fig. 0). 



C. Energy loss 



For many-particle systems, important indicators of the properties of the collision are the net energy and momentum 
losses of the system. Thus we now turn our attention to them. By the time the particles have reached the point of 
rolling without sliding, there has been an energy loss much greater than that which would occur in a system described 
only by a standard coefficient of restitution. For two particles undergoing a collision described by a conventional 
coefficient of restitution, as in Eq. |, the maximum fractional energy loss, which occurs in a head-on collision, is 
(1 — r 2 )/2. For steel balls, with r w 0.903, (1 — r 2 )/2 « 0.09. For collisions that are not head-on, the energy losses 
are smaller, as only the component of velocity normal to the collision decreases, assuming the tangential coefficient 
of restitution equals 1. Figure [s] shows the total system energy versus time for the two-particle collision described 
above. Figure |9| shows the fractional energy remaining in the system at the end of the relaxation time tr for a series 
of collisions as a function of the final angle between the velocities. For a head-on collision the system's energy after 
the collision is on average ~ 37% of the energy before the collision, representing a loss of 63% of system energy as 
the result of a single collision. Also shown in Fig. ^| is a prediction based on Eqs. (solid line) with the parameters 
determined in the fit discussed above. The dashed line represents energy loss in a system with r = 0.903 and no 
surface interactions. The observed energy loss is only weakly dependent on the collision angle for nearly head-on 
collisions. 



D. Momentum loss and the direction of sliding frictional forces 

We note that the direction of force due to sliding friction is not necessarily parallel to the contact normal, h; instead, 
it is in the direction of the relative velocity of the contact between the substrate and the bottom of the particle, which 
we call the contact velocity v ct . Thus, 
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vZt = Vcm + ax LU, (12) 

where v cm is the velocity of the center of mass of the particle, a is the vector from the contact point to the center of 
the particle, and Co is the particle's angular velocity. A sketch is provided in Fig. [l0| Here, we define the x direction as 
Vct/\vct\, and the y direction as (a x i)/|a|. All momentum loss to the substrate will occur in the x direction, as this 
is the direction of the only force acting on the particle (neglecting rolling friction, which is small compared to sliding 
friction). We experimentally determine the x direction by finding the direction in which a particle's velocity changes 
following a collision. Figures |ll| and fl2| show the momenta of the individual particles versus time in the collision 
described above, in the x (Fig. p_l|) and y (Fig. |lj) directions. Note that the x and y directions are independently 
defined for each of the particles, i.e. x for the initially moving particle is different from x for the initially stationary 
particle. Figure |ll| illustrates the finite time after the collision (the collision time t c is marked by a vertical dotted 
line) for which momentum is transferred to the substrate through sliding friction. After this time, the only momentum 
loss is due to rolling friction. In Fig. ^| we observe that no momentum is lost in the y direction for either particle 
after the collision, aside from a slow loss due to rolling friction. 

We also examine the net momentum loss in the direction of the initial momentum versus the final angle between 
the particle velocities (Fig. |l3| ). Note that in the usual case, with no surface interactions, no momentum is lost, so 
p/po = 1 for all angles. In contrast, for a head-on collision with surface interactions we see that ~ 20% of Pi is lost. 
Further, we note that this quantity is weakly dependent on the final angle between the velocities after the collision 
for small angles. The solid line in this figure shows the prediction based on Eqs. |^-^, with the parameters determined 
above. 

This momentum loss may be important in many-particle systems. For example, inelastic collapse, a condition in 
which there are an infinite number of collisions in a finite time, occurs in one- and two-dimensional idealized systems 
p0| , pl| . One-dimensional numerical simulations by Dutt et al. [ ft2[ show that if even a very small momentum loss per 
collision is introduced, inelastic collapse does not occur. This suggests that inelastic collapse cannot be observed in 
experimental granular systems which interact with a surface. 

IV. ANGULAR VELOCITY 

We would also like to determine the angular velocities of the particles immediately after the collision. These are 
difficult to measure directly, but we can derive expressions for them from Eq. [l2[ given the assumption that after 
sliding stops each particle will be rolling without sliding. Then the angular velocity of a particle immediately after 
the collision, uj , is 

Vyo = -(^A« cro + v cmx0 ), (13) 

and 

^xO = --VcmyO, (14) 

where Av cm is the change of the center of mass velocity in the x direction due to sliding forces, v cm o is the center 
of mass velocity immediately following the collision, and v = (1 + ma 2 /I) = 7/2. Since we can directly measure the 
center of mass velocity at all times, we can deduce the values of lu x q and uj v q. 



These calculations determine the components of uj in the x, and y directions, as defined in Sec. [II D. These directions 
are defined by the sliding frictional forces acting on the particles; a more natural coordinate system when examining 
the effect of the collision itself on angular velocities is defined by the n and t directions, that is, parallel (n) and 
perpendicular (t) to the vector connecting the centers of mass of the particles at the time of the collision. If surface 
effects arc negligible during the collision, there is no torque in the n direction, so we expect that Lo n for each particle 
will not be changed by the collision. Indeed, we find that during the collision the mean change of the angular velocity 
in the n direction, averaged over 50 particles and normalized by the angular velocity of the incoming particle in each 
case, is dio n /(v /a) — 0.01 ± 0.02. In contrast, we expect that during the collision, some angular momentum will be 
transferred from the moving particle to the stationary particle in the t direction. The amount of angular velocity 
transferred can be quantified by Eq. [uj we find from these calculations that C = 0.25 ± 0.02. This is similar to, 
although slightly smaller than, the value of C = 0.347 ± 0.008 obtained from the fit to Eqs. above. 
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V. CONCLUSION 



In two-dimensional granular systems, understanding interactions with the substrate is crucial to understanding the 
system dynamics. As two particles collide, there is rotational frustration between them and the substrate, leading 
to sliding on the surface. The large contact force between the particles at the time of the collision is much greater 
than gravity, with (Av / At co u) / g ~ 10 3 , guaranteeing that particles will slide on the substrate after the collision. The 
resulting sliding friction leads to high energy losses and can be modeled simply, as discussed by Kondic In fact, 
we find that up to 63% of the incoming energy is lost in a single collision between two particles with coefficient of 
restitution of 0.9, and most of this is due to sliding friction with the substrate. The sliding continues for a time, 
tr, which is long relative to the collision time, tr is comparable to or longer than the time between collisions for 
moderately dense, rapidly cooling systems, which means that sliding is experimentally important for many-particle 
systems until typical velocities reach v ~ gAt co u . Additionally, we find that both momentum and angular momentum 
are typically lost to the substrate following a collision. 
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FIG. 2. Tracks of two particles colliding: A moving particle enters from the left with v « 10cm/ s and strikes a stationary 
particle. Both exit toward the right. The circles represent the particles' positions every 0.02sec. Inset is a detailed view of the 
particle tracks near the collision point, with the length of one particle diameter shown for scale. 
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FIG. 4. Histogram for the coefficient of restitution, as determined immediately following a collision, for 28 samples. 
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FIG. 5. Angle between particle velocities. The collision occured at t « 0.06 sec, as indicated by the vertical dashed line. At 
this point the particles are moving at nearly 90° to each other. The angle between the incoming particle velocity and h is 5.5°. 
The oscillations are due to experimental noise. 
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□ Particle 1 (initially at rest) 
• Particle 2 (initially moving) 
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FIG. 6. Relaxation time vs. impact velocity. The solid lines are predictions based on a fit to Eqs. |6[j9[ Note that in a system 
without surface interactions, r r = for all collisions. 



11 



35 




q F r i i i 'l' i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i I i i i i 

10 15 20 25 30 35 40 45 50 55 60 

Impact velocity (cm/s) 

FIG. 7. Final velocity vs. initial velocity. The solid lines are predictions based on the fit to Eqs. |^-^| Dotted lines represent 
theoretical final velocities of particles without surface interactions, with r = 0.903. 
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FIG. 8. Total system energy versus time in a single collision. The vertical dotted line marks the approximate time of the 
collision. The dashed line shows the final energy that would result if the fractional energy loss was (1 — r 2 )/2. 
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FIG. 9. Energy loss vs. angle between final velocities of two particles. The dashed line represents numerical calculations of 
two particles colliding without substrate interactions, with r = 0.903. The solid line represents predictions based on Eqs. 
with the fitted parameters r = 0.903, A*fc = 0.232, and C = 0.347. 




FIG. 10. Sketch of velocities for a sliding particle. 
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FIG. 11. Momentum parallel to x vs. time, for (a) the initially stationary particle and (b) the initially moving particle (x is 
different for each particle). The collision occurs at t ~ 0.06sec, as indicated by the vertical dotted lines. Note that a large part 
of the momentum is transfered from one particle to the other at the time of collision, and that v X 2 increases after the collision 
as a result of its spin. 
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FIG. 13. Momentum fraction remaining after relaxation (at time t=t c + T r ) in the direction of the initial momentum (in the 
lab frame) versus final angle between velocities. For a head-on collision, approximately 20% of the total system momentum is 
lost. The solid line gives the prediction based on Eqs. with the fitted parameters r = 0.903, /ife = 0.232, and C = 0.347. 
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